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Abstract 

The pion mass difference generates a pronounced cusp in ^ Svr decays, the 
strength of which is related to the tttt S'-wave scattering lengths. We apply an 
effective field theory framework developed earlier to evaluate the amplitudes for 
Kl — > Stt decays in a systematic manner, where the strictures imposed by analyticity 
and unitarity are respected automatically. The amplitudes for the decay 77 ^ Svr 
are also given. 

Key words: Chiral symmetries, analytic properties of the S'-matrix, decays of 
X-mesons, meson-meson interactions 
PACS: ll.30.Rd, 11.55.Bq, 13.20.Eb, 13.75.Lb 



1. The investigation of the so-called cusp effect in tt+tt^tt^ decays has 

become a fully competitive method for the extraction of the S'-wave vrvr scat- 
tering lengths from experimental data. Following refined versions of the orig- 
inal proposal by Cabibbo [1-3] , the combination oq — 0.2 has been determined 
from very high statistics data [4, 5] to an accuracy mainly limited by remaining 
shortcomings in the theoretical description of the decay amplitudes. Missing 
ingredients are in particular (real and virtual) photon corrections. Here, an 
important step has recently been performed by Isidori [6], who has evaluated 
radiative corrections in multi-body meson decays, in particular, for the fully 
charged channel Tr+Tr+vr", in the soft photon approximation. Once these 

corrections are available in all channels, i^' Svr decays, combined with the 
information gained from K^.^ decays [7, 8] and the pionium lifetime [9] , have 
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the potential to test the very precise theoretical prediction of the scattering 
lengths [10,11] experimentally. For recent phenomenological determinations 
of the scattering lengths, we refer the reader to Refs. [12-14]. 

As the strong impact of the unitarity cusp near the tt'^tt" threshold is a uni- 
versal feature of the n^n'^ scattering amplitude [15], it is present also in other 
decays, like Kl — * 37r°, rj — ^ Svr^ etc. The strength of the cusp in Kl — * Stt^ 
is reduced by about an order of magnitude compared to — > 7t~^7t^t:^, 
hence the experimental situation in order to gain information on tttt scat- 
tering lengths is far less favourable [5] . However, the motivation to study this 
channel all the same is twofold: firstly, experimental efforts to at least see the 
cusp are under way [5]; secondly, the Kl Stt^ system provides an excellent 
object for exploratory studies of the most important electromagnetic effects 
in the cusp region, before immersing oneself into the even more relevant, but 
simultaneously more difficult case of decays. 

The Kl Stt decays have been studied with regard to the cusp phenomenon 
before. Ref. [2] uses unitarity, analyticity and cluster decomposition properties 
of the S'-matrix to investigate the cusp structure. In analogy to the corre- 
sponding K'^ decays discussed in the same reference, an expansion in powers 
of the Tin scattering lengths a is used as the essential ordering principle, and 
the calculation is performed up to O(a^). In Ref. [16], in addition to analytic- 
ity and unitarity, chiral perturbation theory is used for the evaluation of the 
real parts of the K ^ 371 decay amplitudes at one loop. In the present work, 
however, we rely on the non-relativistic effective field theory framework de- 
veloped in Ref. [3]. It is based on an effective Lagrangian, and as such satisfies 
all unitarity and analyticity constraints automatically. The coupling constants 
involved can be directly matched to tttt scattering lengths, and the expansion 
in powers thereof as advocated in Ref. [2] emerges naturally in a generalised 
power counting scheme. 

Our presentation closely follows that of Ref. [3], allowing for a relatively concise 
description of the procedure. We construct the most general non-relativistic 
Lagrangian required for the process in question, and match the couplings to 
the TTTT threshold parameters. Thenceforth the calculation of the decay ampli- 
tude up to two-loop order is straightforward. Our representation of tree, one- 
loop, and two-loop contributions correctly reproduces the analytic structure 
with various branch points and cusps in the Mandelstam plane throughout 
the physical region (and slightly beyond). The pertinent calculation of the 
radiative corrections within the same framework will follow in due course [17]. 

2. We consider the neutral and charged decay modes Kl{Pk) 
7r°(pi)7r°(p2)7r°(p3) and Kl{Pk) {pi)'^ {p2)'^^ {p-i) ■ The kincmatical vari- 

ables are defined as usual: Si = {Pk — PiY with = Mf , i = 1,2,3, 
where Mt^+ = M^^ and M^o denote the masses of the charged and neutral 
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pions, respectively, and A^^ = — M^o 7^ 0. In the centre-of-mass frame 
Pk — {Mk, 0), with Mk the neutral kaon mass, 

2Mk ' ^Ml ' ^ ' 

where A(a;, z) = + + z"^ — 2xy — 2xz — 2yz is the triangle function. 
Below we also use the velocities Vjk and kinetic energies Tj, 

„ , , X(si,Mf,Mi) 
vjkisi)^ ^ f '\ T,^p^-M,. (2) 



3. We invoke the non-relativistic framework set up in Ref. [3] for the evalu- 
ation of the pertinent decay amplitudes. In that framework, the perturbative 
expansion is performed in terms of two formal parameters e and a. One counts 
the pion and kaon masses as 0(1), the pion momenta as 0(e) and the pion 
mass difference A^^ as O(e^). In addition, each four-pion vertex is counted as 
a quantity of order a. As these vertices are proportional to the tttt scattering 
lengths which are small, one expects the expansion in a to converge rapidly. 
We refer for a further discussion of the method to the original article [3]. Here, 
we simply note that it is sufficient to provide the Lagrangian used - the am- 
plitudes then follow from a straightforward application of the rules provided 
in Ref. [3]. 

4. The complete Lagrangian of the effective theory is jCk + ^tt-k: where Cr 
contains Kj^ — > Stt vertices, and vCtttt describes elastic tttt scattering. In the 
following, we provide the Lagrangians necessary to calculate the amplitudes 
for Kl 37r at order e'', ae^, a^e^. 

We start with the tttt interaction and consider the following five physical 
channels in ttV ^ tiH'^: {ab;cd) = (1) (00; 00), (2) (+0;+0), (3) (+-;00), 
(4) (H — ;H — ), (5) (++;++)• [We omit the channel 7r~7r° 7r~7r°, because 
this amphtude is identical to tt+tt^ — > tt+tt^ by charge invariance.] The La- 
grangian takes the form 

^2Y,^iW±{idt - W±)^± + 2^lWo{idt - Wo)^o + E A , (3) 

± i=i 

where <l>j is the non-relativistic pion field operator, W± = ^ — A, Wq — 
^M^o — A, with A the Laplacian. Introducing further the notations 



for n — a, b, c, d, one may write 
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with hi = s\-\ (M^ + M^+M^ + Mj) , where sj denotes the physical threshold 
in the ith channel. Explicitly, hi = 2M^o, /i2 = 2M^M^o, hs = 3M| - M^o, 
hi = = 2M^. The ellipsis stands for terms of order in the S-wave 
and for terms of order in the P- and D-waves. The low-energy constants 
Ci,Di,Ei, Fi are matched to the physical threshold amplitudes below. To sim- 
phfy the resulting expressions, we have furthermore introduced the combinato- 
rial factors xi — — 1/ 4, X2 — X3 — X4 — 1, ui — Us — U5 — 0, U2 — U4 — 1. 
Finally, we note that wc omit local 6-pion couplings. Their contribution to 
the Kl — > Stt amplitude is purely imaginary in the non-relativistic framework, 
and of order e*^. 

5. The couplings C,, A, -E'i, Fi can be expressed in terms of the threshold 
parameters of the underlying relativistic theory. In the isospin symmetry limit, 
the expansion of the relativistic tttt scattering amplitude reads 

+^Af(t -«) + .... (6) 

The eUipsis stands for higher orders in e, e.g. D-wme contributions. The bar 
indicates the isospin symmetric limit, at M^r = 139.57 MeV. In terms of the 
standard scattering lengths qq, 02 and ai, one has 

3^1 = A^(ao + 202) , 2A2^Na2 , 3A3 = A^(a2 - Oq) , 
6^4 = N{2ao + 02) , ^5 = Na2 , 

2A^^Nai , 2A^^Nai , if = if = = , = 327r , (7) 

with ao = 0.220 ± 0.005, = -0.0444 ± 0.0010, Oq - 02 = 0^265 ± 0.004, ai = 
(0.379 ± 0.005) X lO-iM-2 [11]. The products AiU and AJ^ denote effective 
ranges and shape parameters, respectively. 



4 



Still in the isospin symmetry limit, the couplings Ci are related to these thresh- 
old parameters according to 

2Ci = Ai , mlDi = Aifi , 32M^F, = AJi , 4^^ = 3if , (8) 

where we have dropped higher-order terms in the threshold parameters. Tak- 
ing isospin breaking into account, one finds at leading order in chiral pertur- 
bation theory [18] 

2Ci,2,5 = A,2,5(l - V), 2C3 = ^3(1 + r^/3), 2C4 = ^4(1 + V), (9) 

where 77 — ^^./Ml — 6.5 x 10~^. Isospin breaking in the remaining couphngs 
Di,Ei, Fi is expected to have a negligible effect on the analysis, and we propose 
to use for these couplings the relations Eq. (8) also in the real world, where 
isospin is broken. 

6. It remains to display the Kl — > Stt Lagrangian, 
Ck = 2K^WK{idt - Wk)k + Lo(i^^$o^+^- + h.c) 

+Li(k\Wo - M^o) $0*+$- + h.c) + L2(k\Wo - M^of $0*+^- + h.c) 
+L3(i^t$o(W^|*+*- + $+1^1$- - 2W^±$+W^±$_) + h.c) 
+ ^Ko{k^^I + h.c) + ]^K^(k^^I{Wo - M,o)2$o + h.c) + . . . , (10) 
where K denotes the non-relativistic field for the meson, Wk ~ 



Mj^ — A, and the ellipsis stands for the higher-order terms in e. The cou- 
plings Li, Ki are assumed to be real. Their contribution to the decay matrix 
elements at tree-level is provided below. 

The tree-level expressions for the amplitudes, generated by Ck-, are modified 
by final state interactions of the pions, generated by loops evaluated with vCtttt- 
We use the notation 

Mooo = M + M + M 2-l°°P^ + . . . [K:^^ ^W] , 

A^+_0 = A^^^^^ + + + . . . [X^ ^ TT+TT-TT^] (11) 

for the decay amplitudes and the Condon-Shortley phase convention for the 
pions. Our amplitudes are normalised such that the decay rates are given by 

In the case of Kl 37r°, the right hand side must be divided by 3!=6. 
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7. The tree amplitudes are 



7W tree = K, + K, [x^, + Xl + X^j , 
Ml''''' = Lo + + + L3(Xi - X^f , (13) 

where = p° — M^rO . This representation is equivalent to 

■ tree _tt , tt /^„,2 "^^ 



M tree ^ ^ ^ ^^^^^ _ ^ ^^^^^ _ ^ y^^^^ _ ^^y^ ^^4) 
where 

S3 - Sn , f = S2 - Si , 

The relations between the coefficients Ui, Vi and Lj, iCj are displayed in Ap- 
pendix A. 

8. The one-loop contributions are proportional to the basic integral 

"'^ '^J ti2^ 2wa{l)2wb{P - I) K(l)-/o)K(P-l)-Po + /o) ' ^ ^ 

with w±(p) = + p2, -u;o(p) = ^M^o + p2 and ^ p2 _ p2_ ^-^g 

hmit D — >■ 4, 

Ja6(P') = T^^a.(^') , (17) 
iOTT 

which is a quantity of order e. In order to make the formulae more transparent, 
we modify the notation for the couplings Ci, Di, Ei, , 

(Ci, C2, C3, C4, C5) = (Coo, C+o, Cx, C+_, C++) , (18) 

and analogously for the Di,Ei,Fi. In the following, we use J_o = J+o through- 
out, and denote the couplings for 7r~7r° — > 7r~7r° with index -|- as well, C_o = 
C+o, etc. We then find 

^loop ^ !^Bi'\si)Joo{si) + (Si ^ S2) + (Si ^ S3)| 

+ |5j^^(si) J+_(si) + (Si ^ S2) + (Sl ^ S3)|, 

M l-lo«P = 5^(53) Joo (S3) + Bi'Hs3) J+-(S3) 

+ |si')(Si, S2, S3) J+o(si) + (Si ^ S2)|, (19) 
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with 



6si 



Bt^ (S3) = (a + D^Y^, + F^Y^^ \k^ + K 



Xl + 2Zl + ^Y^^ 



B'i\s^) = 2(C+^ + D+^Y,, + F+^Yl) |lo + L1X3 + L^Xl + ^3^>^3c| , 
Si'^ (si, S2, S3) = 2 (C+0 + ^+onnc + i"+onL) + L^Z^ 



+ L2 



(z-)2 + Q? ^23(^1^ 



3si 



+L3 



3si 



1 IT 9l3(si) 



- -E^n^^^l A,(M^ - M|) + s,{s3 - S2) 



x\Li + 2L2Zr + 2L3 Xi -Z+ \ + 0{Al) . 



(20) 



We have used the abbreviations 



Qi-pI+pI (cycl.) , 



2 X{Ml,Mf,s,) 



= - 4M^o , y,e = Si - 4M^ , y,„, = - (m^o + m^)^ , 



(21) 



9. There are two topologically distinct two-loop graphs that describe pion- 
pion rescattering in the final state, see Fig. 1. At the order of accuracy we 
are working, it is sufficient to consider the case of non-derivative couplings. 
In this case, the contributions of both diagrams depend only on the variable 
s, where 

Q'' = (qi + , Q' = s. (22) 

The diagram in Fig. IB, apart from a factor containing coupling constants, is 
given by a product of two one-loop diagrams which were already calculated 
in Eq. (17). The non-trivial contribution from Fig. lA is proportional to 
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Fig. 1. Two topologically distinct non-relativistic two-loop graphs describing the 
final-state tttt rescattering in the decay K — Stt, with = {qi + q2)^- 



(27r)^ i(27r)^ 
1 1 11 



X 



2Wa{\ + k) Wa{\ + k) - Mi^ + /o + fcO 2Wb{\) Wbi}) - P 



111 1 

^ 2wc{k) Wc{k) - P 2wd(Q - k) ^^(Q - k) - go + P ' ^ 



A short discussion of this integral is given in Ref. [3]. There, it is shown that 
one may write 

M{s) = f(m^, M,, M„ M,;s)+..., (24) 



where F is ultraviolet finite and contains the full non-analytic behaviour of 
the two-loop diagram in the low-energy domain, whereas the ellipsis denotes 
terms that amount to a redefinition of the tree-level couplings in and 
which are therefore dropped. A one-dimensional integral representation for F 
is provided in Ref. [3] . The relevant integrals can be performed analytically - 
the result is displayed in Appendix B. 

Below, wc use the notation Fi(. . . ; s) for the integral F{. . . ; s), evaluated at 
Q2 = A(M^, s)/4M^, with i = ±, 0. 

Evaluating the diagrams displayed in Figs. 2 and 3, we find for the amplitudes 
at order a^e^ 

2-loops ^ j^J^A^,^^ + ^B^,^^ ^ ^ ^ ^ ^^)|^ 

2-loops ^ ^A^^^^ ^ ^B^^^^ ^^^^ ^25) 



where 
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Fig. 2. Two-loop graphs contributing to the decay Kl tt vr vr in the non-rela- 
tivistic effective theory. The graphs obtained by a permutation of identical particles 
in the final state are not shown. 



M^(si) = 2 Ko Fo(M^o, M^o, M^o, M^o; si) 
+ 8 C+o a Lo Fo(M^o, M^, M^, M^; si) 
+ 4 Coo Lo Fo{M^, M^, M^o, M^o; si), 



+ 



(2 Cl Ko + 2 Coo a ^o) Joo(5i), 



and 



Miisi, S2, S3) = 2 Coo Ko Fo{M^o, M^o, M^o, M^o; S3) 
+ 4 C^ Lo Fo(M^, M^, M^o, M^o; S3) 
+ 8 C+o C+_ Lo Fo(M^o, M^, M^, M^; S3) 

+ |4 C+o C+_ Lo F+(M,, M,, M^o, M,; si) 

+ 2 C+o C, i^o F+(M,o, M^o, M^o, M,; si) 

+ 4 CIq Lo F+{M^, M^o, M^, M^o; si) + (si ^ S2] 

Mg{su S2, S3) = 4 C^_ Lo 4_(s3) + Coo i^o ^00(^3) 

+ (2 C, C+_ i^o + 2 C^ Lo) J+_(S3) Joo(s3) 

+ |4 Clo Lo Jloisi) + (si ^ S2)|. 



(26) 



(27) 



10. The decay amplitudes depend on the six real Kl — > Stt coupling constants 
Li, Ki and on the threshold parameters for vrvr scattering. Combining the tree- 
and one-loop result Eqs. (13), (19) with the two-loop contributions Eqs. (25), 
we obtain the neutral and charged decay amplitudes up to and including terms 
of order e^, ae^ and a^e^, expressed in terms of the one- and two-loop inte- 
grals J and F displayed in Eqs. (17) and (B.1)-(B.3), respectively. [We have 
dropped some of the contributions at order eA^. In particular, D-waves gen- 
erate contributions of this type. We expect them to be completely negligible.] 
This representation is valid in the whole decay region, and is the main result 
of this article. 
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Fig. 3. Two-loop graphs contributing to the decay Kl 7r"'"7r vr in the non-rela- 
tivistic effective theory. 

The decay amplitude Kl — * n'^n^n^ obeys what we refer to as the threshold 
theorem: the coefficient of the leading non-analytic piece, which is propor- 
tional to ^^-(sa), is given by a product of two factors, the decay amplitude 
Kl —>■ TT^TT+TT" and the scattering amplitude tt+tt" —>■ 7r°7r°, both evaluated 
at threshold [1]. Of course, aside from the determination of the leading term 
in V+-, our approach also allows a systematic evaluation of higher-order con- 
tributions 

11. We now compare the content of this letter with the work of Cabibbo and 
Isidori [2] (CI), who use an alternative method to construct the K ^ 3n decay 
amplitudes. Conceptual aspects of the two methods were already discussed in 
Ref. [3] for the case of the charged kaon decays K'^ Stt. In particular, it 
was pointed out that the amplitudes agree at order a, whereas they differ 
at order away from threshold, because the method used by CI does not 
reproduce the correct analytic properties of the amplitudes at two-loop order. 
[On the other hand, the two amplitudes lead to very similar results for the 
scattering lengths when fitted to K~^ 3n data [5].] Analogous comments 
apply in the case of Kl Stt considered here. Comparing the expressions in 
detail, we note that the final result Eqs. (4.61)-(4.67) in CI does contain some 
(but not all) of the terms evaluated above. In this sense, the expansion of the 
decay amplitudes presented here is more systematic and complete. As to the 
terms retained in CI, we note that, aside from obvious typos, we do agree in 
Kl Stt*^ at order a in the physical region, and at order at the thresholds 
Si = 4M^. In the charged channel Kl — ^ tt+tt^tt^, a graph is omitted in CI. 
It contributes at order ae and generates a cusp at the edge of physical phase 
space. 
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12. We add a remark concerning rj — > 37r° and rj — > 7r"'"7r~7r° decays. These 
processes can be analysed in a completely analogous fashion. Indeed, the tttt 
scattering amplitudes remain the same, whereas the polynomial Lagrangian 
for 77 — > Stt can be obtained from the K — > Stt one by replacing field operators 
and particle masses in the Lagrangian Eq. (10), {K, Mk) {rj, M^). The tree 
amphtudes analogous to Eq. (13) become 

^'^^ = + K1 {Xl + Xl + Xl) , 
Mt'^^ = + LIX^ + VlXl + U>{X, - X,f , (28) 

with Xj = p° — M^o, and with obvious notation otherwise. The relation 
to an alternative expansion in the conventional 77 ^ Stt Dalitz plot vari- 
ables is provided in Appendix A. Furthermore, the one- and two-loop re- 
sults in Eqs. (19), (25) can simply be taken over, with the replacements 
{Ki, Li, Mk) — {Ki,L],M.,j) everywhere. Because r^^^^Tr+Tr-Tro/ri^^^syrO ~ 
r^^7r+7r-7ro/r,j^37rO, wc cxpect that the strength of the cusp effect in the neu- 
tral channel 77 — > 37r° is of the same order as the one in Kl — > Stt", i.e., much 
less visible than in the charged channel K'^ — > tt+tt^tt'^. 

13. In summary, wc have investigated Sn decays within a non- 
relativistic effective Lagrangian framework. The amplitudes are calculated in 
a systematic double expansion in the pion momenta (counted as quantities of 
order e), and in the threshold parameters of elastic tttt scattering (generically 
denoted by a). We provide an explicit representation of the amplitudes at or- 
der e^,ae^,a^e^. The representation is valid in the physical decay region, and 
contains the six (real) — > 3tt coupling constants Lj, and the threshold 
parameters a. The very same amplitude can be used, with trivial modifications 
described above, for a cusp analysis in 77 — > 37r. 

Our amplitudes differ from the ones of Cabibbo and Isidori [2] when compared 
in detail - in particular, we do retain all terms at the above mentioned order 
in the low-energy expansion. For this reason, we believe that it is important 
to check whether our expressions for the amplitudes lead to scattering lengths 
that are in agreement with the ones generated with the amplitudes presented 
in Ref. [2]. 

It remains to investigate radiative corrections, which can be evaluated in the 
field-theoretical framework used here in a standard manner. The effects gener- 
ated by the 7r"^7r~ bound state at the tt'^tt^ threshold can also be investigated 
within the same approach [19-27], see also Ref. [28]. We plan to include these 
effects in forthcoming publications [17]. For the evaluation of radiative correc- 
tions in 7r"'"7r"'"7r~ in the framework of scalar QED, we refer the reader 
to the recent interesting article by Isidori [6] . 
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Appendix A 

The coefficients C/j, Vi are given by 
The inverse relations read 



-s.)\ t/i-^, (A.1) 




K^ = Uo-2U^[{MK-M^of-Sn) , i^i = -Mit/i, 

Lo = Vb + V^[{Mk - M^of - + V2{{Mk - M,.f - s,)' , 

L^=AMkV2{s, - {Mk - A'Uof) - 2MkV, , 

L2 = , L3 = 4MiV3 • (A.2) 



For the decays r] — > 37r° and 77 — > tt'^tt 7r°, the tree amplitudes Eq. (28) may 
be written in the alternative expansion 

M^o ^uo + uiz, Mi*'^^^ ^vo + viy + V2y^ + v^x^ , (A.3) 
with the conventional 77 — > 37r Dalitz plot variables 
^_ v/3(p°-p°) ^ ^_ 3(p°-M.o) ^ ^^2^/3p[-M,y ^ 

where Q?, = — 2M^ — MttO, Qj^o = -^?? ~ "iM^a. The coefficients of the two 
representations are related by 
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vi = ^ + Iq,LI) , V, = , V, = ^Ll , (A.5) 



or reversely by 



K^^Uo- 2ui , Kl = -pr^Ui , L2^Vo-Vi+V2 , 



6 

QriO 



Appendix B 

The analytic expression for the two-loop function F reads 

F{M,,M,,M„Ma,s) = N {2A h + B f^) + 0{e^) , (B.l) 



with 



/o = 4(t;i + ^2 - ^2 + , 



/l = ^ - 1) + y2{v2 - 1) - y2(^^2 - 1) + /i) , 

Vi — V~yi arctan —= , i — 1,2 ; V2 — V~y2 arctan ■ 



-B T VB^ - AAC 

yi,2 = ^ , ^2 = ^/2 [S) 

s s 
.0 = Mi + M'^ - 2Mk{mI + Q!il±i)!j ' ' , 

^,_ \{Ml,Ml{M, + M,f) ^_ Ml-Ml 

4M|- ' s ' ^ ' ' 
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The arctan is understood to be evaluated according to 

1 1 + ix /T-. \ 

arctanx = — m , (B.3) 

and s is given a small positive imaginary part in all arguments, s — > s + ie. 

The analytic formula Eq. (B.l) is exact at O(e^), and thus at the order con- 
sidered in Ref . [3] . It differs by a few percent from the integral representation 
given in Ref. [3]. 
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